calculate phase transition